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Abstract 

The lithium problem arises from the significant discrepancy between the primordial 7 Li abun- 
dance as predicted by BBN theory and the WMAP baryon density, and the pre-Galactic lithium 
abundance inferred from observations of metal-poor (Population II) stars. This problem has loomed 
for the past decade, with a persistent discrepancy of a factor of 2 — 3 in 7 Li/H. Recent developments 
have sharpened all aspects of the Li problem. Namely: (1) BBN theory predictions have sharpened 
due to new nuclear data, particularly the uncertainty on 3 He(a,7) 7 Be, has reduced to 7.4%, and 
with a central value shift of ~ +0.04 keV barn. (2) The WMAP 5-year data now yields a cosmic 
baryon density with an uncertainty reduced to 2.7%. (3) Observations of metal-poor stars have 
tested for systematic effects, and have reaped new lithium isotopic data. With these, we now find 
that the BBN+WMAP predicts 7 Li/H = (5.24tg;^) x 10~ 10 . The Li problem remains and indeed 
is exacerbated; the discrepancy is now a factor 2.4 - 4.3 or 4.2a (from globular cluster stars) to 
5.3<7 (from halo field stars). Possible resolutions to the lithium problem are briefly reviewed, and 
key nuclear, particle, and astronomical measurements highlighted. 
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I. INTRODUCTION 



Big bang nucleosynthesis (BBN) remains one of few probes of the early Universe with 
direct experimental or observational consequences lfl. Historically, BBN was generally taken 
to be a three-parameter theory with results depending on the baryon density of the Universe, 
the neutron mean-life, and number of neutrino flavors. Indeed, concordance between theory 
and observation for the abundances of the light elements, D, 3 He, 4 He, and 7 Li was a 
powerful tool for obtaining the baryon density of the Universe. Over the last twenty or so 
years, the number of light neutrino flavors has been fixed (in the Standard Model of particle 
physics) j^j], and the neutron mean-life is measured accurately enough that it is no longer a 
parameter, but rather its residual uncertainty is simply carried into primarily an uncertainty 

n 

in the predicted He abundance (and is small) [2J. 

More recently, the baryon density has been determined to unprecedented accuracy by 
WMAP pj to be Vl B h 2 = 0.02273 ± 0.00062 where £l B = p B /p c is the fraction of critical 
density in baryons, p c = 1.88 x 10~ 29 h 2 g cm -3 , and h is the Hubble parameter scaled 
to 100 km/Mpc/s. This corresponds to a baryon-photon ratio of rj w = 6.23 ± 0.17 where 
V — %/w 7 = 10~ 10 ?7io. This is significantly more accurate than any determination of rj from 
observational determinations of light element abundances. 

Thus the patadi gm for BBN has shifted in the post-WMAP eta fl Q. Nue.eosyttthesis 
is now a parameter-free theory. In the Standard Model, and working in the framework of a 
Friedmann-Robertson- Walker cosmology, the only inputs are the nuclear reaction rates (and 
their associated uncertainties). The last significant update to many of the needed rates was 
compiled by the NACRE g roup 6] and recent BBN calculations by several groups are in 



good agreement 



lid 



ii 



12|. 



What has emerged is rather excellent agreement between the predicted abundance of D /H 
as compared with the determined abundance from quasar absorption systems 13] . Indeed 
what is often termed the success in cosmology between BBN and the CMB is in reality 
only the concordance between theory and observation for D/H at the WMAP value of rj. 
Currently, there is no discrepancy between theory and observation for 4 He. But this success 
is tempered by the fact that 4 He is a poor baryometer, it varies only logarithmically with rj, 



and the observational uncertainty it rather large 



14 



15|. 



It has also become generally accepted that there is a problem concerning the abundance 
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of 7 Li. At the WMAP value of i], the predicted abundance of 7 Li is approximately 3 times 
the observationally determined value. Several attempts at explaining this discrepancy by 
adjusting some of the key nuclear rates proved unsuccessful {^[lf], 17]. In fact, the key pro- 
cess for the production of 7 Li (in the form of 7 Be) is 3 He(a, 7) 7 Be and has been remeasured 
recently by several groups 18j and one of the purposes of this paper is to explore the conse- 
quences of these new rates on the BBN production of 7 Li. At this time, it is unclear whether 
or not any known systematic effect may be responsible for the discrepancy or whether there 
is an indication of new physics at play. 

In what follows, we will review the status of the 3 He(a,7) 7 Be reaction rate and employ 
a new thermal fit to this rate. We will also briefly review the status of the 7 Li observations 
as well as that of the other light elements. In section III, we will present the results of 
the BBN calculations using the new 3 He(a,7) 7 Be rates and quantify the discrepancy with 
observations. In section IV, we discuss the remaining known alternative explanations of the 
discrepancy and we give our conclusions in section V. 



II. NEW LIGHT ON THE LITHIUM PROBLEM 
A. BBN Theory: Updated Nuclear Data 

The possibility of systematic errors in the 3 He(a,7) 7 Be reaction, which is the most im- 
portant 7 Li production process in BBN, was considered in detail in [17|. The absolute value 
of the cross section for this reaction is known relatively poorly both experimentally and 
theoretically. However, the agreement between the standard solar model and solar neutrino 
data provides additional constraints on possible systematic shifts in this cross section. Us- 



ing the standard solar model of Bahcall 19J, and recent solar neutrino data 20j, one can 
exclude systematic variations of the magnitude needed to resolve the BBN 7 Li problem at 



the > 95% confidence level 



In order to maintain up-to-date nuclear input for BBN, several new efforts need to be 
addressed. In particular, the recent experiments measuring the neutron lifetime, evaluations 
of the p(n, 7)<i and 3 He(a,7) 7 Be data, and the resulting revisions to the 3 He(a,7) 7 Be cross 
section and its error budget. 

n 

The new neutron lifetime measurement by |2l|] finds a lifetime, r n = 878.5 ± 0.8. This is 
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more than 5<r away from the current world average r n = 885.7 ± 0.8 recommended by the 
Particle Data Group [2[. Obviously, this points to unknown systematics between the exper- 
iments. Several possible explanations exist 22j], but remain to be explored experimentally. 
We thus, follow the PDG'08 recommendation, until some understanding of the discrepancy 
can be found. 

The p(n, 7)<i determines at what point the deuterium bottleneck ends. Experimentally, 
this reaction is quite difficult to measure in the relevant energy range for BBN. We thus rely 
on theoretical-based cross-sections which are normalized to the existing experimental data. 



Several methods 
potential models 



rave been used to determine the shape of this cross section; microscopic 



23J], i?-matrix fits [2J], N 4 LO pion-less effective field theory [2a], and NLO 



di-baryon effective field theory 



261 ] . All methods agree well, though the R- matrix fit deviates 



by up to ~ 4% from the others. Evaluating all the relevant data, np-scattering, deuteron 

ground state properties as well as the radiative capture and photo-dissociation data using 

a Markov Chain Monte Carlo algorithm, a best fit with propagated uncertainties 26( , finds 

< 1% errors in the capture rate. We adopt this rate in our network. The new rate is ~1% 

□ 

higher than that adopted by [10J at E ~ 500 keV. 
At 

^wmapi 7 Li is produced as 7 Be, and the only important reactions are production via 
3 He(a,7) 7 Be and destruction via 7 Be(n,p) 7 Li and subsequent 7 Li(p, a) 4 He. Of these, new 
developments in 3 He(a, 7) 7 Be experimental efforts [l8j, must be taken into account. The ap- 
parent discrepancy between prompt gamma-ray and 7 Be activation measurements discussed 
in [27j has virtually disappeared, though discrepancies between datasets still persist. A new 
evaluation 28], examines the "modern" data and finds a best fit and 68.3% CL uncertainties 
for a model independent fit, using a Markov Chain Monte Carlo algorithm similar to that 



of 



261 ] . Special care is taken to properly propagate the dominating systematic normalization 



errors, as standard statistical treatments underestimate the true uncertainties. They find 
a zero-energy astrophysical S factor S^O) = 0.580 ± 0.043 keV b, an uncertainty of 7.4%. 
In the Gamow window (E ~ 200 keV), this new fit gives a central value 17% higher than 
previous work 101 ] . 

n 

Using the scalings provided in [10J (their Eqn. 47), the two new rates and baryon density, 
we expect an increase in the 7 Li prediction to ~ 5.3 x 10 -10 . This increase combined with 
reduced errors will aggravate the Li problem as we show below. 
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B. CMB Cosmological Parameters: Updated Baryon Density 



As noted above, the major paradigm shift for BBN came with the first WMAP de- 
termination of the baryon density. The first year WMAP best fit, assuming a vari- 



able spectral index was 29] fish 2 = 0.0224 ± 0.0009 which corresponds to a value of r\: 
7/io(WMAP2003) = 6.14 ± 0.25. This level of accuracy was already able to pin the light 
element abundances down to a narrow range. 

The 5-year data WMAP data are consistent with their first year data, and the errors 
have been significantly reduced. The 5-year data give fish 2 = 0.02273 ± 0.00062, or 

7? 10 (WMAP2008) = 6.23 ± 0.17 (1) 

We will adopt this value for the present study. 

It is worth recalling key assumptions which enter into the inference of rj from the CMB. 
The baryon density arises from the pattern of acoustic oscillations which modulate an un- 
derlying primordial spectrum of temperature perturbations (e..g., [2^]). The baryon density 
inferred from CMB data thus depends on the nature of the primordial perturbation power 
spectrum, which is usually assumed (with motivation from inflation) to be a simple or run- 
ning power law, close to the Harrison-Zel'dovich form. If the primordial power spectrum 
had a richer form, however, then the cosmological parameters inferred from the tempera- 
ture anisotropies can be changed substantially 3l|. Fortunately, the inclusion of precision 
polarization data, and their cross-correlation with the temperature data, goes far to in- 
dependently characterize the primordial spectrum and to sharpen the precision of cosmic 
parameters. The launch of Planck in the next few months will offer just such precision deter- 
mination of CMB polarization, and thus should further reduce the systematic and statistical 
uncertainties in the baryon density. 

In using the WMAP value for rj at the period of BBN, we are implicitly assuming that 
there has been no entropy change between BBN and the decoupling of the CMB. Note that 
entropy production between BBN and decoupling would require a larger value for 77 at the 
time of BBN and make the Li problem even worse. 
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C. Stellar Lithium Observations: Updates, Systematics and Isotopic Ratios 



There have been several recent observational determinations of the 7 Li abundance in 
metal-poor halo stars. Most observations lead to a 7 Li abundance in the range (1 — 2) x 10 -10 , 



consistent with the original determination by Spite and Spite 



suggested a small but significant correlation between Li and Fe 33] which can be understood 



3211 . Precision data have 



as the result of Li production from Galactic cosmic rays {34. I35 ] . Extrapolating to zero 
metallicity one arrives at a primordial value |36|] Li/H| p = (1.23 ±0 06) x 10" 10 , though the 
systematic uncertainties were recognized to be large and an abundance of 



(Li/H) 



p, field* 



;i.23±£g) x 10- 



10 



(2) 



was derived at 95% confidence. 

There have been several measurements of 7 Li abundances in globular clusters. The abun- 
dance in NGC 6397 was measured by to be 



(Li/H) P)G c = (2.19 ± 0.28) x 10 



-10 



(3) 



Other groups found similar values of (1.91 ± 0.44) x 10~ 10 [23] and 7 Li/H = (1.69 ± 0.27) x 



10 10 39j. A related study (also of globular cluster stars) gives 7 Li/H = (2.29 ± 0.94) x 
10" 10 4(J- 

An important source for systematic error lies in the derived effective temperature of the 
star. [Li] = log( 7 Li/H) + 12 is very sensitive to the temperature, with d[Li]/dT c g ~ 0.065 
- 0.08. Unfortunately there is no industry standard for determining effective temperatures, 
and for a given star, there is considerable range depending on the method used. This 
spread in temperatures was made manifest in the recent work of Melendez and Ramirez 41 ] 
using the infra-red flux method (IRFM) which showed differences for very low metallicities 
([Fe/H] < -3) by as much as 500 K, with typical differences of ~ 200 K with respect to 
that of [3^ 1 . As a consequence the derived 7 Li abundance was significantly higher with 
Li/H| p = (2.34 ±0.32) x 10 



-10 



41 



42 1 



Recent observations of metal-poor stars do not find support for the large temperature 
scale in {4]]]. Using H-a wings Asplund et al. |4^] found temperatures consistent with 
previous determinations using the IRFM. This work also showed a non-negligible correlation 
with metallicity leading to an estimate of the primordial 7 Li abundance between 1.1 and 



7 



1.5 xlO~ 10 . A larger slope (with larger uncertainty) was also seen in [44] (also using H- 
a lines) in their study of extremely metal-poor stars. Because of the large slope seen in 
this data, the extrapolation to zero metallicity leads to a primordial 7 Li abundance of 8.72 
±1.71 x 10~ n , below that found in previous studies (though their mean value is consistent 
with other determinations). 

A dedicated set of observations were performed with the specific goal of determining 
the effective temperature in metal-poor stars [2] . Using a large set of Fe I excitation lines 
(~ 100 lines per star), the Boltzmann equation was used with the excitation energies, Xi 
to determine the temperature through the distribution of excited levels. Again, there was 
no evidence for the high temperatures reported in 41], rather, temperatures were found 
to be consistent with previous determinations. The mean 7 Li abundance found in [45j was 
Li/H = (1.3 — 1.4±0.2) x 10~ 10 , consistent with the bulk of prior abundance determinations. 

There are of course other possible sources of systematic uncertainty in the 7 Li abundance. 
It is possible that some of the surface 7 Li has been depleted if the outer layers of the stars 
have been transported deep enough into the interior, and/or mixed with material from the 



hot interior; this may occur due to convection, rotationa 
possible depletion factors are in the range ~ 0.2-0.4 dex 



mixing, or diffusion. Estimates for 



461 ] . Recent attempts to deplete the 



Li abundance through diffusion introduce a source of turbulence tuned to fit the abundances 



of heavy elements in NGC6397 



471 ] . It is not clear whether this mechanism will work for 



the wide range of stellar parameters seen in the field. As noted above, the Li data show a 
negligible intrinsic spread in Li. Any mechanism which reduces significantly the abundance 
of Li must do so uniformly over a wide range of stellar parameters (temperature, surface 
gravity, metallicity, rotational velocity etc.). At the same time, the mechanism must avoid 
the nuclear burning of 7 Li in order to preserve some 6 Li (which burns at a lower temperature). 

In fact, not only does 6 Li challenge models of stellar depletion of 7 Li, but may req uire non- 
standard models to understand its abundance at low metallicity. Recent data 43J (though 
called into question [48[]) indicates the presence of a plateau in 6 Li at the level of about 
1000 times the BBN predicted 6 Li abundance 35|, |49|]. While this abundance of 6 Li can 



not be explained by conventional galactic cosmic-ray nucleosynthesis 



34!, 



35 



501 ]. it can be 



explained by cosmological cosmic-ray nucleosynthesis due to cosmic rays produced at the 



epoch of structure formation 



5l|. 



It is also possible that the lithium discrepancy(ies) is(are) a sign of new physics beyond 
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the Standard Model. One possibility is the cosmological variation of the fine structure 
constant. Varying a would induce a variation in the deuterium binding energy and could 
yield a decrease in the predicted abundance of 7 Li 52j . A potential solution to both lithium 
problems is par ticle decay after BBN which could lower the 7 Li abundance and produce some 
6 Li as well [53J. This has been investigated in the framework of the constrained minimal 
supersymmetric Standard Model if the lightest supersymmetric particle is assumed to be 



;he gravitino 



55]. 



54| and indeed, some models have been found which accomplish these goals 



III. THE LITHIUM PROBLEM QUANTIFIED 

As discussed in section IIA, we have incorporated the new cross section measurements 
of 3 He(a,7) 7 Be as well as an updated rate for p(n,7)d into the BBN calculation of the 
light element abundance. Since the Li problem is our primary focus here, we first show the 



effect of the new rate determination for 3 He(a,7) 7 Be 28| in Fig. [Q There, the primordial 
abundance of 7 Li/H is plotted as a function of r\ or equivalently the baryon density, Q b h 2 , 
shown on the upper horizontal axis. The new result is shown by the green shaded region 
which is superimposed on top of the older result taken from [10|] and shown shaded in red. 
The thickness of the bands corresponds to one sigma uncertainty in the calculated abundance 
as determined by a Monte Carlo simulation of the BBN reaction network. Also shown as a 
vertical yellow shaded band, is the one sigma range for 77 as determined by WMAP [3]. To 
highlight the change in the Li error budget, Fig. [2] shows the fractional uncertainty in Li/H 
as a function of 77, for the same two sets of rates as in Fig. [TJ 

As expected, at low rj, there is virtually no effect of the new rate for 3 He(a,7) 7 Be, with 
essentially no change in the central value and fractional error of 7 Li. However, at higher 77, 
Fig. [1] shows a slight shift to larger values of 7 Li/H. Most pronounced at high 77, however, is 
the reduction the 7 Li fractional error seen in Fig. [2J indicating a greatly reduced uncertainty 
in the BBN calculation. At the WMAP value of r] 10 = 6.23, we find 

7 Li/H = (5.24ijS) x 10- 10 (4) 

to be compared with the previous value of 7 Li/H = (4.26toJg) x 10 -10 [10|]. The shifts in 
both the central value and the error range almost entirely reflect the corresponding shifts in 
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baryon density Vl b h z 
0.01 




10 



-9 



baryon-to-photon ratio 77 



FIG. 1: A comparison of the 1A/H abundances as a function of n using the new 3 He(^a, 7,) 7 Be rate 



lu l (red band). The thickness of the bands 



determined in J 28] (green band) to the older result from 
represents 1 a uncertainties in the calculated abundance. Central values are given by thick solid 
and dashed curves respectively. The thin upper curve (within the green shaded region) demarcates 



the border of the old result. The yellow vertical band is the WMAP value 



ofvm- 



the 3 He(a,7) 7 Be cross section and rate. 

One may notice the difference between the baryon-density convolved value of 
7 Li/H=4.26 x 1(T 10 (rj 10 = 6.14 ±0.25) and the fixed baryon density value of 7 Li/H=4.36 x 
10~ 10 (7710 = 6.14). This ~ 2% shift is due to 77 dependence of the uncertainties in the 7 Li/H 
prediction. With the higher precision year 5 WMAP baryon density, this shift is < 1%. The 
differences can be fully explained by the simple scalings from jlo| . 

In {9], it was found that the adopted prescription for defining errorbars, underestimated 
the true uncertainty for the reaction 3 He(a,7) 7 Be, due to the visible discrepancy between 
datasets. Therefore, a different method was used to evaluate this reaction's uncertainty (and 
only this reaction). Had this method been uniformly applied to all reactions, to account for 



10 



all discrepant or non-discrepant data, recognizing that the quoted errors are in fact minimum 
(not actual) uncertainties akin to the "minimal" errors discussed in [7J, their prediction would 
still agree with the result in [10| . but with inflated errorbars ( 7 Li/H — 4.15^^5 x 1(T 1U ). 



baryon density Cl h h z 
0.01 




baryon-to-photon ratio 77 

FIG. 2: The fractional error <r(Li)/Li in the BBN lithium prediction, shown as a function of n. 
The solid curve is our result; the broken curve gives the older result from \lC J as in Fig. [TJ The 
yellow vertical band is the WMAP value of r\ 

The full set of light element abundances from BBN is shown in Fig. [3j The 4 He abundance 
is plotted as a mass fraction as a function of 77 in the top panel. The thickness of the 4 He band 
is primarily due to the small uncertainty in the neutron mean life. The D and 3 He abundances 
by number with respect to H are shown in the middle panel, and the 7 Li abundance (also 
by number) is shown in the lower panel. As one can see, the WMAP determination of the 
baryon density is precise enough, that one can now simply read off the predicted primordial 
abundances of the light elements. These are: 

(5) 
(6) 



Y p = 0.2486 ± 0.0002 



D/H = (2.49 ±0.17) x 10" 



11 



and 



3 He/H = (1.00 ±0.07) x 10~ 5 

baryon density n.h 2 

y 10-2 y b 



(7) 




FIG. 3: The light element abundances of D, 3 He ; 7 Li by number with respect to H, and the mass 
fraction o/ 4 He as a function of r\. The thickness of the bands represents 1 a uncertainties in the 
calculated abundance. The yellow band gives the WMAP r\ Jj/. 



The BBN predictions can be compared directly with current observational determina- 
tions of the light element abundances. The BBN likelihood functions can be defined by a 
convolution over r\ 

^bbnPO = dr) L BB n(ri\X) L WMAP (r)) (8) 
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using the Monte Carlo results from BBN as a function of r\ to give Lbbn(?7|^0 an d the 
WMAP value of r\ distributed as a Gaussian, Lwmap(^)- These are shown in Fig. H]by the 
dark (blue) shaded regions. Though there are useful measurements of the 3 He abundance 



56(, these are difficult to match to the primordial abundance 57|. We will show the BBN 



likelihood for 3 He in Fig. HI but will not discuss 3 He any further. 
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FIG. 4: The theoretical and observational likelihood functions for 4 He, D/H, 3 He/H, and Li/H. 
BBN results have been convolved with the WMAP determination of rj and are shown as dark (blue) 
shaded area. The observational likelihoods are shown as light (yellow) shaded regions as well as 
alternative dotted curves. The data and distinctions are detailed in the text. 



Fig. H] also shows the observational likelihoods for comparison. For 4 He, the light (yellow) 
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shaded region corresponds to the result found in [14| using a select subset of the data in 



58(. Systematic uncertainties and degeneracies in the set of physical parameters used to 



determine the helium abundance led to a value of Y extrapolated to zero metallicity of 
Y p = 0.249 ± 0.009. The mean value of the reanalyzed data is 0.252 ± 0.003 and that is 
shown by the dotted likelihood function. Because of the large observational uncertainty in 
determining Y p , both agree quite well with the accurate prediction made by BBN. 

The deuterium abundance at low metallicity has been measured in several quasar absorp- 
tion systems 13] . The weighted mean value of the seven systems with reliable abundance 
determinations is log D/H = —4.55 ± 0.03 where the error includes a scale factor of 1.72 [2| 
and corresponds to D/H = (2.82 ±0.21) x 10~ 5 . Since the D/H shows considerable scatter, 
it may be appropriate to derive the uncertainty using sample variance (see eg. 7|) which 
gives log D/H = -4.55 ± 0.08 or D/H = (2.82 ± 0.53) x 10~ 5 . Both are shown plotted as 
log-normal distributions by the dashed curve (light shaded area) and dotted curve respec- 
tively in the D/H panel of Fig. HI As one can see, there is very good agreement between the 
predicted and observationally determined value of D/H. 

Finally, we come to 7 Li. The inferred primordial value from field halo stars (eq. [2]) is 
shown by the light shaded area in the Li panel of Fig. HI The slightly higher determination 
in a globular cluster (eq. [3]) is shown by the dotted likelihood function. As one can see, 
the reduced uncertainty in the BBN prediction leads to virtually no overlap between the 
theoretical and observational likelihoods. 



IV. EXPERIMENTAL AND OBSERVATIONAL STRATEGIES 

We see that the cosmic lithium problem not only remains but has become even more 
pressing. As we have noted already, several proposed solutions to the Li discrepancy remain 
viable. Fortunately, the different scenarios proposed thus far can be ruled in (or out) empir- 
ically, both in the laboratory and in the observatory. Here we outline a strategy for doing 
so. 

The most exciting possibility, in our opinion, is that the lithium problem points to new 
physics at work in the early Universe. This possibility could receive strong support-in the 
most optimistic case, nearly outright confirmation-by the discovery of new physics at the 
Large Hadronic Collider (LHC) at CERN, slated to start running in late 2008. The presence 
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of particle physics beyond the Standard Model would immediately cast the Li problem in a 
new light, since the correctness of the Standard Model is one of the key assumptions in the 
BBN calculation of Li and all of the light element abundances. 

Consider, for example, the implications of the LHC finding evidence for supersymmetry 
(SUSY). This would place the highest priority on the study of SUSY effects on nucleosyn- 
thesis based on the potentially measurable sparticle spectrum. However, the fact that super- 
symmetry is realized in nature would not, by itself, guarantee a solution to the Li problem. 
This also requires that SUSY particles alter the light element abundances either during or 
after BBN. The lightest supersymmetric particle (LSP) would be stable and make up the 
weakly-interacting dark matter today, but its feeble interactions would also make it an in- 
ert spectator to nucleosynthesis. If, however, the next-to-lightest supersymmetric (NLSP) 
particle were relatively long lived (r ~ 1 — 1000 sec or more) then its decays would release 
a cascade of Standard Model particles which would alter the light element abundances and 



potentially solve the Li problem 53|, |54|, |55j . Hence, the most direct laboratory signature of 



a SUSY solution to the Li problem would be the discovery of a long-lived NLSP. 

In many of the cases studied 54|, |55[, the gravitino is the LSP and the partner of the tau 
lepton is the NLSP. Unfortunately, in the constrained versions of the minimal supersym- 
metric Standard Model, if these two sparticles make up the low mass end of the spectrum, 
obtaining lifetimes as short as ~ 1000 s, would require a very massive spectrum beyond the 
reach of the LHC. 

New astronomical observations of various kinds will play an important role in further 
sharpening the Li problem and probing its solutions. To date, metal-poor halo stars in 
our own Galaxy remain the only observational targets for determining the primordial Li 
abundance. Identifying and exploiting new sites for primitive Li detection, with different 
systematic issues, would offer crucial cross-checks on the halo star results. For example, 
metal-poor high velocity clouds probably represent the Galactic infall of the primitive intra- 
group medium of the Local Group. Li abundances in these clouds, seen in absorption against 
background quasars, would provide the first measures of extragalactic Li and with very deep 
exposures could even yield isotopic abundances, all without any issues of stellar evolution 
or atmospheric modeling 5j|. An even more challenging but exciting possibility would be 



to observe highly redshifted lines from cosmic Li recombination 



60J. This process occurs 



after the usual hydrogen recombination but still at z ~ 400 — 500 and thus the Li features, 
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if measured, would probe not only the cosmic Li abundance before star formation began, 
but also give information about the universe in the "dark ages." Another observation which 
bears indirectly on the Li problem will be the new determination of the extragalactic gamma- 
ray background by the now-operational GLAST observatory; these data probe the non- 
cosmological but pre-Galactic Li contribution by cosmic rays in galaxies and on cosmological 



scales 



If and when new sites for Li determination become available, halo stars will nevertheless 
remain the dominant probe of primordial Li for some time to come. Thus there will remain a 
pressing need to understand systematics and to obtain abundance measures which are both 
precise and accurate. Further targeted, systematic studies of the kind outlined in §11 CI are 
of the highest importance. Moreover, robust 6 Li data are exceedingly important, as the mere 
existence of primordial 6 Li at currently observable levels immediately demands nonstandard 
nucleosynthesis, and the 6 Li/H abundance is a crucial constraint on all models. Even in the 
absence of primordial 6 Li, the inevitable component from Galactic cosmic-ray processes, and 
a possible contribution from cosmological cosmic rays, both provide unique and empirical 
means of subtracting the non-primordial 7 Li produced by these mechanisms. 

Nuclear experiment and theory remain important inputs to the Li problem, but as noted 
above, it would seem that the few relevant reactions are by now well understood, and a 
nuclear solution to the Li problem is unlikely. That said, it is obvious that the discovery of a 
large and unknown systematic error would not only be surprising but could also dramatically 
affect the problem. 

Finally, we note that these diverse experimental and observational arenas will contribute 
in complementary ways. As LHC probes of SUSY and other new physics rule particle 
solutions in or out, astronomical observations will sharpen the ability of Li and BBN to 
probe the solution space. 



V. CONCLUSIONS 



The first-year WMAP data determined the cosmic baryon-to-photon ratio r\ with a high 
precision that, in concert with BBN theory, predicts a primordial 7 Li abundance significantly 
above the levels inferred from the most metal-poor halo stars. In this paper we have revis- 
ited this cosmic lithium problem, which has seen important new contributions from nuclear 
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experiments, halo star observations, and the 5-year WMAP data. We find that the discrep- 
ancy is now a factor 2.4- 4.3 or 4.2a (from globular cluster stars) to 5.3a (from halo field 
stars), which is largely due to the tighter errors on the baryon density and the 3 He(a, 7) 7 Be 
cross section. 

A wide variety of possible solutions to the Li problem have been proposed. These range 
from observational systematics in determining halo star Li, to nuclear effects, to exotic 
physics beyond the Standard Model. Of these, nuclear physics solutions seem essentially 
ruled out at this point. Observational systematics remain a challenging problem and a pos- 
sible solution, but the hope that a revised temperature scale in halo stars would significantly 
reduce the Li discrepancy is challenged by recent data which confirm the existing scale. The 
possibility that lithium points to new physics at work in the early universe thus remains not 
only viable, but if anything more likely. 

Upcoming experiments and observations will help clarify the lithium problem and its 
solution. In particular, the LHC will probe particle physics Li solutions generally and su- 
persymmetry in particular. These results will be complemented by a host of astronomical 
observations. Planck observations of the CMB, particularly polarization, will further nail 
down the cosmic baryon density. The first extragalactic Li abundances will offer an inde- 
pendent constraint on primordial Li with completely different systematics than halo stars. 
And finally further halo star measurements, focusing on systematics and on 6 Li, will remain 
crucial in quantifying the Li discrepancy from standard BBN and probing the mechanism 
which resolves this disagreement. 

R.H.C is supported through the NSF grant PHY 02 16783 (JINA). The work of KAO 
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